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Abstract. We discuss the potential of a next generation space-borne CMB experiment for 
studies of extragalactic sources with reference to COrE+, a project submitted to ESA in 
response to the call for a Medium-size mission (M4). We consider three possible options 
for the telescope size: lm, 1.5 m and 2 m (although the last option is probably impractical, 
given the M4 boundary conditions). The proposed instrument will be far more sensitive 
than Planck and will have a diffraction-limited angular resolution. These properties imply 
that even the 1 m telescope option will perform substantially better than Planck for studies 
of extragalactic sources. The source detection limits as a function of frequency have been 
estimated by means of realistic simulations taking into account all the relevant foregrounds. 
Predictions for the various classes of extragalactic sources are based on up-to-date models. 
The most significant improvements over Planck results are presented for each option. COrE-f- 
will provide much larger samples of truly local star-forming galaxies (by about a factor of 8 
for the lm telescope, of 17 for 1.5m, of 30 for 2m), making possible analyses of the prop¬ 
erties of galaxies (luminosity functions, dust mass functions, star formation rate functions, 
dust temperature distributions, etc.) across the Hubble sequence. Even more interestingly, 
COrE-f- will detect, at \b\ > 30°, thousands of strongly gravitationally lensed galaxies (about 
2,000, 6,000 and 13,000 for the lm, 1.5m and 2m options, respectively). Such large samples 
are of extraordinary astrophysical and cosmological value in many fields. Moreover, COrE+ 
high frequency maps will be optimally suited to pick up proto-clusters of dusty galaxies, i.e. 
to investigate the evolution of large scale structure at larger redshifts than can be reached 
by other means. Thanks to its high sensitivity COrE+ will also yield a spectacular advance 
in the blind detection of extragalactic sources in polarization: we expect that it will detect 
up to a factor of 40 (lm option) or of 160 (1.5m option) more radio sources than can be 
detected by Planck and, for the first time, from several tens (1 m option) to a few hundreds 
(1.5 m option) of star forming galaxies. This will open a new window for studies of the global 
properties of magnetic fields in star forming galaxies and of their relationships with SFRs. 

Keywords: cosmology: observations - surveys - submillimeter: galaxies - radio continuum: 
general - galaxies: evolution 

1 Introduction 

We investigate the impact on studies of extragalactic sources of the planned fourth gener¬ 
ation Cosmic Microwave Background (CMB) mission named Cosmic Origins Explorer plus 
(COrE+). Various options are being considered. The COrE+ “light” concept envisages an 
effective telescope size of ~ 1 m and a 60-600 GHz frequency range with a total of 2040 
detectors. More ambitious options, requiring substantive contributions from international 
partners (“COrE+ extended”), contemplate telescope sizes from ~ 1.5m (baseline) to 2m 
(although the last option is unlikely to fit within the M4 constraints) and a frequency range 
from 60 to 800 GHz (but we consider also the effect of an extension to 1200 GHz) with a 
total of 5800 detectors. The instrument comprises more than two times as many frequency 
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Figure 1. Left panel: power spectra of the astrophysical components present in COrE+ maps at 
70 GHz at high Galactic latitudes (|6| > 30°), computed from the Planck Sky Model [2], masking the 
region at |6| < 30° with a 10° apodization for 30° < |6| < 40°. The three dot-dashed straight green 
lines show the power spectra of radio sources, that are the dominant extragalactic population at this 
frequency, computed after having masked sources brighter than the detection limits corresponding 
to the three options for the COrE+ telescope, shown in Fig. 4. The contributions to the power 
spectrum of Galactic emission are detailed in the adjacent panel. Right panel: power spectra of 
the astrophysical components present in COrE+ maps at 150 and 353 GHz. The lines have the same 
meaning as in the left panel. However at these frequencies the only relevant Galactic emission is the 
one of thermal dust. Also, fluctuations from dusty galaxies, which have a large contribution from 
clustering, are increasingly important with increasing frequency. 


channels as Planck , to make possible an efficient separation of foreground emissions. The 
very large number of state-of-the-art detectors ensures a much better sensitivity per channel. 

The plan of the paper is the following. In Section 2 we discuss the source detection in 
CMB maps. In Section 3 we briefly describe the populations of extragalactic sources in the 
C0rE+ frequency range. Section 4 is about the extraction of the rich information content on 
sources below the detection limit. In Section 5 we deal with counts in polarization. Finally, 
in Sect 6 we summarize our main results. 

Throughout this paper we adopt the fiducial ACDM cosmology with best-fit Planck + 
WP + high-resolution CMB data of the Atacama Cosmology Telescope (ACT) and South 
Pole Telescope (SPT) experiments as provided by [1], 


2 Detection limits for CMB maps 
2.1 Basics 

The development of algorithms for point source detection has a long history in astronomy. 
Images produced by CMB experiments have from this point of view important peculiarities 
[3] that make inadequate the algorithms consecrated by decades of successful usage in many 
wavebands, from radio to X-rays, such as CLEAN [4], DAOFIND [5] and SExtractor [6]. The 
peculiarities mostly follow from the limited angular resolution of CMB experiments. 

In high resolution images intensity peaks on the angular scale of the Full Width at 
Half-Maximum (FWHM) of the telescope are due to either point sources or to detector noise 
and the highest significance peaks are straightforwardly interpreted as point sources. On the 
contrary, in CMB images also structure in diffuse emissions (the CMB itself and Galactic 
radiations) shows up on the FWHM scale, complicating the detection of point sources. 
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Figure 2. Contributions of the various astrophysical components to the root mean square (rms) 
intensity fluctuations at high Galactic latitude (|6| > 30°) at the angular resolutions of the three 
options for the COrE+ telescope, as a function of the observed frequency. Fluctuations are dominated 
by the CMB up to 271 GHz and by Galactic thermal dust at higher frequencies. The dominance of 
diffuse emissions highlights the importance of efficient component separation algorithms. The WMAP 
and Planck experiences have demonstrated that fluctuations due to diffuse emissions can be largely 
filtered out by source extraction algorithms. The contributions of extragalactic sources have been 
computed after masking sources brighter than the detection limits shown in Fig. 4. 


Optimal methods for source detection in CMB maps need the knowledge of the proba¬ 
bility density functions (PDFs) of each signal contributing to the “background” [for a review 
see 3]. In general, however, our knowledge of the PDFs is highly incomplete. A further com¬ 
plication comes from the non-stationarity of Galactic emissions, hampering the realization 
of detection methods valid throughout the sky. Still, methods able to efficiently filter out 
the “background” noise have been devised [7]. The widely used matched filter characterizes 
the “background” noise in terms of its power spectrum. Although this approach has known 
limitations [3] it constitutes a useful reference. 

The power spectra at high Galactic latitude (\b\ > 30°) of the astrophysical components 
present in COrE+ maps at selected frequencies (70, 150 and 353 GHz) are illustrated in 
Fig. 1. They were computed using the most recent version of the Planck Sky Model [2] 
with an apodized Galactic mask (see the caption of Fig. 1). These figures illustrate the 
relative importance of the various contributions to the fluctuation field and how they vary 
with frequency. 

The root mean square (rms) fluctuations, cr, are related to the power spectrum, C £ , by 
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is the solid angle of the instrument, b £ is the experimental beam function, which for a Gaussian 
beam with standard deviation 9 b is well approximated by 


b £ = exp[-^0^ + 1)]. 
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The integral approximation to the sum holds if 9 b <C lrad (so that the relevant l 1). 
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Figure 3. Root mean square (rms) intensity fluctuations at the resolution of the measurements as 
a function of the full width at half maximum (FWHM) of the telescope. The deviation of Poisson 
fluctuations from the direct proportionality, a oc FWHM, suggested by the combination of eqs. (2.6) 
and (2.11), is due to the decrease of the detection limit, Sd, in eq. (2.6) as the angular resolution 
improves. The decrease of the fluctuation amplitude with decreasing Sd is stronger for flatter slopes 
of the source counts. 


Poisson fluctuations are uncorrelated. Hence the only non-zero term of the correla¬ 
tion function is the one at zero-lag and the power spectrum, which is the transform of the 
correlation function, is independent of £: 
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where dN(S)/dS are the differential number counts per steradian of sources weaker that the 
detection limit Sd- Then 
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which is the classical expression for the variance of a Poisson distribution of sources weaker 
than Sd, within a solid angle cj e ff defined by 


CU e ff 27T 



d eefiOA). 


(2.7) 


Adopting a Gaussian instrumental response function (the “beam” function), in polar 

coordinates: 
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Figure 4. PCCS 90% completeness limits [thick red 
bars; 9] compared with 1 a fluctuations due to CMB 
and Galactic emissions. The latter have been esti¬ 
mated using the mask described in the caption of 
Fig. 1. The fact that the detection limits are gen¬ 
erally not far from such rms fluctuations, while the 
threshold signal to noise ratios for source detection 
range from 4 to 4.9, illustrates the power of source ex¬ 
traction algorithms in filtering out diffuse emissions. 
The estimated 4 a detection limits for the 3 COrE-l- 
options (dot-dashed green lines) are also shown for 
comparison. The dashed orange curve shows the 
5cr noise levels for the Planck nominal mission (15 
months). 
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where we have used the relation FWHM = 2\/2 In 2 0^ ~ 2.3550;,. Note that the solid angle to 
be used to compute Poisson fluctuations is half of that to be used in the general case (eqs. 2.1 
and 2.2). 

The dependence of rms fluctuations on frequency for diffraction limited observations 
with a 1 m, 1.5 m and a 2 m telescope is illustrated by Fig. 2. Figure 3 shows the trend of rms 
fluctuations with angular resolution at 3 frequencies. The relative importance of fluctuations 
due to unresolved point sources increases with frequency and with decreasing FWHM. Their 
dominant contribution in the FWHM range of interest here comes from radio sources up to 
about 150 GHz and from dusty galaxies at higher frequencies. While the radio sources have 
an essentially Poisson distribution, in the case of dusty galaxies the contribution of clustering 
dominates except on the smallest angular scales. 

The above discussion illustrates why, in spite of the very high sensitivity of modern 
CMB experiments, they only provide shallow surveys of extragalactic sources. For example, 
the 90% completeness limits of the Planck Catalogue of Compact Sources [PCCS; 9] at 70, 
100, 217 and 857 GHz are 776, 300, 180 and 680 mJy, respectively, while the 5cr instrumental 
noise limits at the same frequencies, for the 15 month survey, are ~ 235, 131, 143, and 220 
mJy, respectively [10, 11]. This means that the surveys are limited by fluctuations of sky 
signals, not by instrumental noise: there is a lot of useful information in the fluctuation held! 
A comparison between the PCCS 90% completeness limits at all Planck frequencies and the 
corresponding 5a instrumental noise levels is presented in Fig. 4. 

A further cue comes from the consideration of the number of detections per resolution 
element. For example, the PCCS has 191, 629, 1409 and 6773 detections at 70, 100, 217 
and 857 GHz, respectively, above the 90% completeness level in the extragalactic zone. This 
means that there are 1940, 1119, 1868 and 516 resolution elements per source, well above the 
values (30-40 resolution elements/source) usually corresponding to the confusion limit. 

A key issue is the ability of source extraction algorithms to filter out fluctuations due 
to diffuse emissions, such as the CMB and Galactic radiations. Figure 4 shows that the 
component separation algorithms developed for the Planck project [7, 12] are quite efficient 
at doing that: over a broad frequency range the PCCS 90% completeness limits are close to 
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the 1 a fluctuations due to diffuse emissions. In other words most of their contributions to 
fluctuations are filtered out. 

2.2 Estimate of COrE+ detection limits 

To estimate the detection limits for the three COrE+ options considered in this paper we have 
carried out realistic simulations of sky maps including all the relevant foregrounds as well 
as the instrumental noise for a three year mission. Simulations with resolutions > 4arcmin 
have been carried out using the latest version of the Planck Sky Model [PSM: 2]. The PSM 
includes, in addition to a model of the CMB, Galactic diffuse emissions (synchrotron, free- 
free, thermal and spinning dust, CO lines), Galactic HII regions, extragalactic radio sources, 
dusty galaxies, thermal and kinetic Sunyaev-Zeldovich signals from clusters of galaxies. Each 
component is simulated by means of educated interpolations/extrapolations of the available 
data, complemented by state-of-the-art models. 

Distinctive features of the simulations are spatially varying spectral properties of syn¬ 
chrotron and dust; different spectral parameters for each point source; modelling of the 
clustering properties of extragalactic sources and of the power spectrum of fluctuations in 
the cosmic infrared background (CIB), in close agreement with the latest observational de¬ 
terminations. 

On angular scales below 4arcmin, reached at frequencies above ~ 160-320 GHz (de¬ 
pending on the telescope size), the PSM has not been well tested yet. So, at high frequencies 
we have performed simpler simulations, taking advantage of the fact that, on small scales, 
the fluctuation field at high Galactic latitudes is dominated by extragalactic sources making 
up the CIB and fluctuations due to the CMB and to diffuse Galactic emissions can safely 
be neglected. The method developed by [13] was used to distribute the sources on the sky 
consistent with the measured CIB power spectra at different sub-mm wavelengths [14-16], 
i.e. properly allowing for their clustering properties. 

The source detection on simulated maps was performed using the IFCAMEX detec¬ 
tion pipeline 1 and specifically its implementation of the second member of the Mexican Hat 
Wavelet family (MHW2). This wavelet is obtained applying the Laplacian operator two 
times on the two-dimensional Gaussian function. The MHW2 is used as a filter to reduce the 
contribution from the background, including both the small scale noise and the large scale 
diffuse emission from our own Galaxy, thus enhancing the detection efficiency. The analysis 
was carried out in projected flat 3.66° x 3.66° square patches, corresponding 128 x 128 pixels 
for HEALPix 2 [Hierarchical Equal Area isoLatitude Pixelation; 17] maps with A r s id e = 2048. 

For this purpose the map was divided into partially overlapping flat patches. Dividing 
the sky map into small patches we can optimize the filter taking into account the statistical 
properties of the map in the vicinity of each source. In practice, we did first a blind run 
and then a second non-blind run at the position of each source detected in the blind run. 
The second step allowed us to remove artifacts in the blind catalogue, mostly due to border 
effects in the filtered images. 

Since the flat patches overlap, we can have multiple detections of the same source. The 
repetitions are removed, keeping for each source only the detection with the highest signal- 
to-noise ratio. Further details on the MHW2 and the IFCAMEX implementations can be 
found in [18] and in [19]. The IFCAMEX detection pipeline was previously used to build the 
New Extragalactic WMAP Point Source (NEWPS) catalogue [20, 21] and the part of the 

1 http: //max.ifca.unican.es/IFCAMEX 

2 http://healpix. sourceforge.net 
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Figure 5. Spectral energy distributions (SEDs) of a 
bright radio source ([HB89] 0234+285; dashed black 
line) and of a bright infrared galaxy (M61; solid black 
line) compared with the estimated 4 a detection lim¬ 
its for the three COrE+ telescope sizes (dot-dashed 
green lines). The radio emission of M61 was es¬ 
timated exploiting the mean relationships between 
IR and radio emissions derived by [26]. The syn¬ 
chrotron and free-free contributions are shown by 
the dashed black lines at the bottom of the fig¬ 
ure. Synchrotron dominates at low frequencies. Note 
that the frequency at which the radio spectrum of 
[HB89] 0234+285 crosses the dust emission of M61 
is only 0.3 dex higher than the crossing frequency of 
its own radio emission, which is almost 2 orders of 
magnitude fainter. 


PCCS at v < 70 GHz; it also being used to build the second Planck catalogue of compact 
sources (PCCS2) at the same frequencies. A different implementation of the same algorithm 
was used for the PCCS and PCSS2 catalogues at higher frequencies. 

The 4cr detection limits, Sd, that, according to our simulations, correspond to > 90% 
completeness and reliability in the “extragalactic zone” (|6| > 30°), where Galactic emissions 
are quite low, are shown in Figs. 4 and 5. The values of Sd obtained from the simulations 
are well approximated by the formula: 


Sd — 4[(7c 0n f + 0-noise + (0.12 ffCMB)"]^, 


( 2 . 10 ) 


with 
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In the above formulae the various contributions to the fluctuation field have been computed 
using the equations given in Sect. 2.1. The CMB power spectrum corresponds to the Planck 
best fit cosmological parameters [1]. The rms confusion fluctuations, cr con f, include the con¬ 
tributions of radio sources, of dusty galaxies and of galaxy clusters (Sunyaev-Zeldovich, SZ, 
effect). The Poisson contribution from radio sources was computed via eq. (2.6) using the 
[22] model at frequencies of up to 100 GHz and the [23] model at higher frequencies. The 
detection limit, Sd, was computed iteratively. A check with the PSM has shown that the 
clustering of radio sources adds a negligible contribution (we recall that the PSM contains, 
down to faint flux density levels, real radio sources, at their real positions in the sky). The 
contributions of dusty galaxies and of galaxy clusters were computed using, respectively, the 
model by [24], that fits accurately the measured CMB power spectra, and the power spectrum 
of the thermal SZ effect measured by [25]. These power spectra comprise both the Poisson 
and the clustering contributions. 


3 Extragalactic sources in the COrE+ frequency range 

By a lucky (for CMB studies) coincidence the CMB peak occurs at a frequency close to 
a deep minimum of the intensity of the extragalactic background light as well as of the 
Galactic emissions. The steep increase with frequency of the dust emission spectrum in the 
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Figure 6. Differential number counts of radio sources at 70, 150, 220 and 550 GHz. The Planck data 
points are from [27] at 70 GHz and from [29] at higher frequencies (counts based on the PCCS are 
not available yet; the 90% completeness limits of the latter catalogue in the “extragalactic zone” are 
indicated by the vertical solid lines). The South Pole Telescope (SPT) points are from [28] and the 
Herschel point at 550 GHz is from [30]. The curves are the model predictions by [22] at 70 GHz and 
by [23] at higher frequencies. The vertical dashed lines correspond to the 4 a detection limits for the 
3 COrE+ telescope options. 



Figure 7. Differential number counts of dusty galaxies at 150, 220, 550, 850 and 1200 GHz (2 mm, 
1.4 mm, 545 fim, 353 /im and 250 pm). The Planck ERCSC data points are from [29] and [31] . The 
SPT points are from [28] and the Herschel points from [32]. The curves are the model predictions 
by [24]: the red lines show the contributions of unlensed (solid) and strongly lensed (dashed) proto- 
spheroidal galaxies, the cyan lines show the counts of late-type normal and starburst galaxies, the 
back lines are the total. The vertical dashed lines correspond to the 4 a detection limits for the 3 
COrE+ telescope options while the vertical solid lines indicate the 90% completeness limits of the 
PCCS in the “extragalactic zone”. 


mm/sub-mm region (typically S v oc is 3 ' 5 ) while the radio emissions (synchrotron and free- 
free) decline makes the crossover frequency between radio and dust emission components 
only weakly dependent on their relative intensities. Moreover, dust temperatures tend to be 
higher for distant high luminosity sources, partially compensating for the effect of redshift. 
As a consequence there is an abrupt change in the populations of bright sources above and 
below ~ 1 mm: radio sources dominate at longer wavelengths, while in the sub-mm region 
dusty galaxies take over (Fig. 5). 

The Planck Early Release Compact Source Catalogue [ERCSC; 33] reported sources 
extracted from maps built using the data obtained from the scans of the sky between 2009 
August 13 and 2010 June 6, comprising the first all-sky survey and about 60% of the second. 
It was exploited by [27] to derive counts of extragalactic radio sources in the range 30-217 
GHz and by [29] to investigate the statistical properties of infrared and radio extragalactic 
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sources between 100 and 857 GHz. 

Subsequent releases of the Planck compact source catalogue benefited of additional data 
(15 months, “the nominal mission”, for the PCCS, full mission for the PCCS2) as well as of 
improved data processing. As a result, the completeness limits of the PCCS are fainter than 
those of the ERCSC by average factors of ~ 1.5 at LFI frequencies (30-70 GHz) and of ~ 2.5 
in the range 100-353 GHz; the improvement is larger at higher frequencies: a factor of 3.5 
at 545 GHz and of 4.5 at 857 GHz. No information is available yet on the PCCS2. However 
major improvements are not expected. 

Differential counts at several frequencies of radio sources and of dusty galaxies are shown 
in Figs. 6 and 7, respectively. The Planck counts (yellow data points) were obtained from 
the ERCSC (the PCCS has not been exploited yet for this purpose). The 90% completeness 
limits of the PCCS and of the 3 COrE+ options are also shown. The scientific advances 
made possible by COrE+ are discussed in the following sub-sections. 

3.1 Classical radio sources 

COrE+ covers the frequency range where the information on the spectral energy distribution 
(SED) of radio-loud Active Galactic Nuclei (AGNs) is still scanty: even the SEDs of very 
bright sources like [HB89] 0234+285 (Fig. 5) are not measured at sub-mm wavelengths. Yet, 
in this spectral region important spectral features, carrying essential information on physical 
conditions of sources, show up. 

Observations at mm/sub-mm wavelengths often reveal the transition from optically 
thick to optically thin radio emission in the most compact regions, i.e. on the maximum self¬ 
absorption frequency which roughly corresponds to the synchrotron peak frequency (in terms 
of vL u ). A systematic survey in the COrE+ range will, for example, allow us to see if there 
are systematic differences in the synchrotron turnover frequencies between BL Lacs and flat- 
spectrum radio quasars, as would be expected if their jets have different distributions of the 
angular separations from the line of sight, implying different amounts of Doppler boosting. 
Correlations between turnover frequency and luminosity, which is also boosted by relativistic 
beaming effects, would help confirm current models. 

Major high radio frequency flares have been observed in several compact radio sources 
[34-36], including the recent giant outburst of 3C 454.3 [37]. Several mechanisms can yield 
strong variability of relativistically beamed sources [blazars; e.g. 35, and references therein]: 
shocks travelling along the jet, changes of the bulk Lorentz factor causing variations of 
the Doppler factor, changes of the viewing angle due, e.g., to precession of binary black 
hole systems or to helical trajectories of plasma elements or to rotating helical jets. Thus 
variability studies provide clues into the physical properties of emitting regions. Variability 
due to geometrical effects is expected to affect weakly if at all the source spectra, while 
shocks induce strong spectral variations with intensity peaks that generally move downwards 
in frequency with shock age. The COrE+ surveys may catch the rise of the flare at the 
highest frequencies, missed by ground based observations. 

As for extended sources, the spectral break frequency, i/&, at which the synchrotron 
spectrum steepens due to electron energy losses, is related to the magnetic field, B, and 
to the synchrotron age, t s (in Myr), by Vb — 96(30/iG/R) 3 t~ 2 GHz. The systematic multi¬ 
frequency study at the COrE+ frequencies will thus provide a statistical estimate of the 
radio source ages. On the other hand, at high redshifts the dominant energy loss mechanism 
of relativistic electrons may be inverse Compton scattering off CMB photons because the 
CMB energy density, which grows as (1 + z) 4 , can exceed the magnetic energy density in 
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Blazar SEDs 



Figure 8. SEDs of three blazars detected by 
Planck showing a sub-mm excess attributable to 
host galaxies endowed with active star formation. 
Since, with the exception of galaxies with extreme 
SFRs, Planck could only detect the far-IR emis¬ 
sion from nearby galaxies, the sub-mm detection 
of these blazars implies that a substantial fraction 
of these sources are hosted by star-forming galax¬ 
ies, contrary to the notion that there host galaxies 
are passive ellipticals. The Planck data are repre¬ 
sented by the black points with error bars; the other 
data points were taken from the NASA/IPAC Ex- 
tragalactic Database, NED. For clarity, the SEDs of 
GB6 1347+1217 and of NGC 1275 have been raised 
by A log {vf v ) = 3 and = 6, respectively. 


radio lobes. As a consequence the synchrotron emission of extended, steep-spectrum, sources 
is suppressed with increasing redshift down to lower and lower frequencies, while compact 
sources are almost unaffected [38] . This has obvious implications for the statistics of extended 
sources as a function of redshift and of the survey frequency. Detailed predictions depend on 
poorly known quantities, such as the distribution of magnetic field intensities, of source sizes 
and of injected relativistic electrons. A comparison of low-frequency ground based surveys 
with the high frequency surveys by COrE+ can provide key constraints on these quantities. 

Excess far-IR/sub-mm emission due to dust mostly heated by young stars is often ob¬ 
served from radio galaxies, particularly at high z [39]. COrE+’s broad frequency coverage 
will allow an extensive study of the spectral energy distribution from radio to sub-mm wave¬ 
lengths and therefore of the presence of dust in host galaxies and of the relationships between 
the radio emission and the evolutionary status of the host galaxies. 

According to the current notion, host galaxies of blazars are passive ellipticals [40]. 
This view has however already been challenged by Planck observations of catalogued blazars 
associated to galaxies with intense star formation activity (3 examples are shown in Fig. 8). 
Since dust emission from the host galaxy is detectable by Planck only for nearby galaxies 
(but blazars are rare locally) or for extreme IR luminosities, these objects are likely only the 
tip of the iceberg. The higher sensitivity of COrE+ will allow us to extend the study to 
much larger samples, shedding light on the possible relation of the blazar phenomenon with 
star-formation activity. 

3.2 Special radio source populations 

High radio frequency surveys are crucial to investigate special classes of radio sources, self- 
absorbed up to cm wavelengths. A particularly important class of sources in this category are 
extreme GHz Peaked Spectrum (GPS) sources or high frequency peakers [HFPs; 41]. GPS 
sources are powerful [log(Li. 4 GHz/ergs -1 Hz” 1 ) £ 32], compact (sizes ;$ lkpc) radio sources 
with convex spectra peaking at GHz frequencies. Conclusive evidence that they correspond 
to the earliest stages of the evolution of powerful radio sources, when the radio-emitting 
region grows and expands within the interstellar medium of the host galaxy came from VLBI 
measurements of propagation velocities of up to ~ 0.4c, implying dynamical ages of ~ 10 3 
years [see 42, for a review and references]. The identification of these sources is, therefore, a 
key element in the study of the birth of radio activity. 
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Predictions for the counts of GPS sources have been worked out by [22, 43] and [45]. 
Observational estimates are complicated by the contamination of candidate GPS samples by 
flaring blazars that also exhibit spectral peaks at GHz frequencies [44]. Recent estimates 
differ by substantial factors, but the uncertainties are very large due to the poor statistics. 
[46] found that only 3 sources with S 20 GHZ = 200 mJy selected over an area of about 6.1 sr 
could be classified as genuine GPSs, with a possible incompleteness by a factor of 2. On the 
other hand, [47] have 2 confirmed GPSs with S 20 GHZ > 200 mJy, both with u p > 80 GHz, over 
a much smaller area; this would correspond to 24-155 (68% confidence limits for a Poisson 
distribution) GPSs over 6.1 sr. Only all sky surveys can produce significant samples of these 
rare objects, but they need to go deeper than Planck. COrE+, especially with the largest 
telescope size, has the right properties. 

3.3 Local dusty galaxies 

The (sub-)mm surveys performed by all-sky CMB experiments are the ideal way to carry out 
an unbiased census of dusty galaxies in the local neighbourhood, down to volume densities 
beyond the reach of more sensitive surveys over small areas of the sky, such as those carried 
out by Herschel. The Planck ERCSC has offered the first opportunity to accurately determine 
the luminosity functions in the very local Universe at several (sub-)millimetre wavelengths, 
unaffected by cosmological evolution [31], and to investigate their properties [48]. These 
studies have not yet been extended to the PCCS. 

In Sect. 2 we have shown that, at sub-mm wavelengths, COrE+ will reach, at 500- 
600 GHz flux densities a factor of ~ 4, 6.7 and 9.6, for the lm, 1.5 m and 2 m option, 
respectively, fainter than the PCCS, i.e. to explore, respectively, a volume a factor of 8, 17 
and 30 times larger. This means that COrE+ can detect at 600 GHz, from ~ 7, 900 (1 m 
option) to ~ 17,800 (1.5 m option) to ~ 31,500 (2 m option) star forming galaxies out to 
0 = 0.1 in the “extragalactic zone” . 

This redshift range is within that covered by the Sloan Digital Sky Survey (SDSS) whose 
DR10 release extends over 14,555 square degrees [49]. Additional redshifts are provided by 
other, albeit shallower, large-area spectroscopic redshift surveys: the all-sky 2MASS Redshift 
Survey [2MRS; 50] and the IRAS PSC Redshift Survey [PSCz; 51]; the hemispherical Six- 
Degree Field Galaxy Survey [6dFGS; 52], and others, such as the Two-Degree Field Galaxy 
Survey [2dFGRS; 53]; see, e.g., [54] for a compilation of redshift measurements, dubbed 
2M++. 

Accurate photometric redshifts have been obtained for many more galaxies detected by 
wide-angle surveys [55-57]. The 2MASS Photometric Redshift (2MPZ) sample of [55] should 
be particularly useful in this respect due to its almost full sky coverage, suitable depth and 
size (1 million galaxies with a median z = 0.07; for comparison, the median values of the 
low -0 peaks of the distributions in Fig. 9 are ~ 0.02-0.03). In addition, COrE+ local dusty 
galaxies missed by 2MASS should be present in the WISE catalogue [58]; a large fraction, if 
not all of them, will have reliable photometric redshifts available over most of the sky [59, 
Bilicki et al. in prep.]. We expect that with the new data that will be accumulating in the 
coming years, particularly those from the Euclid slitless spectroscopy and from 4-metre Multi- 
Object Spectroscopic Telescope (4MOST) multifibre R = 5000 resolution spectroscopy of 10 
million z < 1 emission line galaxies, spectroscopic or photometric redshifts will be available 
for essentially all the 0 < 0.1 galaxies detected by COrE+. 

The large number of galaxies for which spectroscopic or photometric redshifts will be 
available will make possible analyses of the properties of galaxies (luminosity functions, dust 
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Figure 9. Predicted redshift distributions of dusty galaxies at 295, 390 e 600 GHz at the 4 a detection 
limits of the 1 m, 1.5 m and 2 m options, from top to bottom. Such limits are 81.0, 44.1 and 29.1 mJy 
at 295 GHz; 104.5, 60.7 and 41.6 mJy at 390 GHz; 150.5, 90.4 and 63.5 mJy at 600 GHz. The distri¬ 
butions are strikingly bi-modal with two well separated peaks. The low-z peak is due to relatively 
local (z 0.1) star-forming galaxies while the high-z (z £ 1.5) peak is dominated by strongly lensed 
galaxies. 


mass functions, star formation rate functions, dust temperature distributions, etc.) across 
the Hubble sequence. Of special interest will be the study of IR emission from IR-faint 
galaxies such as early-type or dwarf galaxies. Combining the COrE+ with the available 
multi-wavelength data, especially with those from WISE, 2MASS, GALEX and Euclid we 
will get a complete view of the SEDs and will be able to investigate relationships between 
the dust content, the star formation rate (SFR), the stellar mass, the environment and more. 

For dusty galaxies, that dominate the extragalactic source population at high frequen¬ 
cies, lower frequency measurements provide information on the relationship, if any, between 
star-formation and nuclear radio activity and on the radio emission (synchrotron and free- 
free) powered by star formation. 
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3.4 High redshift dusty galaxies 

Sub-mm surveys [see, e.g., 62, for a recent review] have clearly demonstrated the importance 
of the FIR/sub-mm data in reconstructing a complete picture of the history of galaxy for¬ 
mation and evolution. The high-redshift sources detected in these surveys are expected to 
be the progenitors of the present-day massive elliptical galaxies [63-65]. These galaxies can 
reach extreme IR luminosities [> lO 13 Z/ 0 ; 66-69], that can be further boosted by strong 
gravitational lensing [70-74] to the point that some were detected by Planck, in spite of the 
shallowness of its surveys [75-79]. 

Using the [24] model, that accurately fits a broad variety of infrared to millimeter- 
wave data on extragalactic sources 3 (multi-frequency and multi-epoch luminosity functions 
of galaxies and AGNs, redshift distributions, number counts, total and per redshift bins) we 
find that the deeper COrE+ surveys, compared to Planck, will detect unlensed IR galaxies up 
to redshifts > 1 (for telescope sizes > 1.5 m) and many more strongly lensed (amplification 
/i > 2) galaxies. In the “extragalactic zone” (\b\ > 30°) we expect, at 600 GHz, ~ 540 
unlensed galaxies at z > 1 for the 1.5 m option. For the 2 m option this number jumps to 
~ 15,900. 

The counts of strongly lensed galaxies were computed using again the [24] model, but 
adopting a maximum amplification /r max = 30 to take into account the sources size [at 
variance with 24, who assumed point sources]. As shown by Fig. 2 of [80] the adopted /i max 
is consistent with the counts and the redshift distribution of strongly lensed galaxies detected 
by the South Pole Telescope [SPT; 28, 81]. 

Even the 1 m option will detect thousands of strongly lensed galaxies at |b| > 30° (about 
740 at 390 GHz and over 2,000 at 600 GHz). For the 1.5m option these numbers increase 
to ~ 2, 900 and ~ 6 , 200, respectively; for the 2 m option they are ~ 6 , 700 and ~ 12, 700, 
respectively. 

For all the options considered the COrE+ counts will cover the gap between the bright 
(sub-)mm counts measured by Planck and those measured by the deep SPT and Plerschel 
surveys (orange and green points, respectively, in Figs. 6 and 7). The transition from the 
bright counts, dominated by nearby galaxies, which have a Euclidean slope, to the fainter 
ones made much steeper by the combined effect of strong cosmological evolution and of the 
positive K-correction [82, 83] provides a strong test for galaxy evolution models. Many of 
them predict a smoother transition than suggested by present data, but only the much better 
statistics that will be provided by COrE+ will make it possible to reach firm conclusions. 

The large samples of strongly gravitationally lensed galaxies that will be provided by 
COrE+ will be of extraordinary astrophysical value in many fields [see 84, for a review] 
thanks to the magnification of the source flux that makes observable galaxies intrinsically 
too faint, and to the corresponding increase of the apparent size that makes possible to 
measure the internal structure of high-z sources to levels otherwise unattainable with the 
current instrumentation [e.g., 85]. Follow-up observation will allow us to determine the total 
(visible and dark) mass of the lensing galaxy, to investigate its density profile, to measure 
cosmological parameters [86-89] and especially the Hubble constant using time delays [e.g., 
90]. 

The selection of strongly lensed galaxies will be extremely easy at about 300-400 GHz. 
As illustrated by Fig. 9, essentially all high-z galaxies detected at these frequencies will 
be strongly lensed. The other detected sources will be easily recognizable low-z late-type 

3 See figures in littp://staff.ustc.edu.cn/~zcai/ 
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Figure 10. Distribution of 
intensity peaks of the Cosmic 
Infrared Background (CIB) 
at 545 GHz at the Planck res¬ 
olution (5'), computed using 
the PSM. A Gaussian fitting 
its negative side is also shown 
for comparison. 
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galaxies, plus a small group of radio sources, also easily identifiable in low frequency radio 
catalogs. At 600 GHz the fraction of strongly lensed galaxies among z > 1 objects is still very 
close to 100% for the lm option and only slightly decreases (to 92%) for 1.5 m option. For 
the 2 m option the strongly lensed fraction among 600 GHz detections with z > 1 is 44%; but 
in this case we expect ~ 6, 700 strongly lensed galaxies in the “extragalactic zone” already 
at 390 GHz, where they are ~ 87% of objects at z > 1. This means that COrE+ will allow 
us to find very easily thousands of strongly lensed galaxies with a close to 100% efficiency 
[see also 72], 

While other facilities will also be generating large reliable gravitational lens catalogues 
on a comparable timescale [e.g. Euclid Ha lenses [91], Gaia lensed quasars [e.g. 92]], the 
critical advantages of COrE+ will be in extending the sources and lenses to much higher 
redshifts. This will make it possible to probe the evolution of dark matter halo substructure 
and the Initial Mass Function (IMF) to much higher redshifts [key tests of semi-analytic 
models of galaxy evolution, e.g. 93, 94]. This large, high-redshift lens catalogue will also be 
ideal to probe cosmological parameters using rare double lenses [see, e.g., [95] for the use 
of dynamical constraints to remove model degeneracies, e.g. [96]]. The background source 
population will mostly be ultra-luminous star-forming systems that in themselves have always 
posed the strongest challenges to semi-analytic galaxy evolution models. Only in large lens 
samples can the rarest high-magnification events be found, and the angular magnification 
afforded by gravitational lensing in these systems will make them the ideal laboratories for 
determining the physical processes that dominate star formation and feedback in the early 
Universe. 

4 Extracting information on extragalactic sources below the detection 
limit 

4.1 Proto-clusters of dusty galaxies 

As mentioned in Sect. 2, thanks to the very low instrumental noise, the fluctuation field 
measured by COrE+ is signal dominated. This implies that it contains a lot of useful infor- 
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mation. In regions with low dust content, several intensity peaks not associated to individual 
sources just below the detection limit were found to correspond to clumps of dusty galaxies 
likely to evolve into rich galaxy clusters [76, 79, 97]. This was predicted by [98] on the basis 
of the argument summarized and slightly updated below. 

Sub-mm surveys proved to be most efficient in detecting high -2 massive galaxies, in¬ 
terpreted as progenitors of present day giant ellipticals, caught during their star-formation 
phase. There is evidence of strong clustering of these sources [14-16, 99-101], consistent with 
them being tracers of strongly overdense regions. Sub-mm surveys are therefore optimally 
suited to look for proto-clusters at earlier redshifts than can be reached by optical/near-IR, 
X-ray, SZ surveys. 

The analysis in Sect. 2 has demonstrated that, in high Galactic latitude regions, intensity 
peaks below the point source detection limits at sub-mm wavelengths are dominated by 
fluctuations in the distribution of faint sources making up the CIB. Such fluctuation field is 
highly skewed. One example, computed using the PSM at the Planck resolution, is displayed 
in Fig. 10. The clustering contribution overcomes Poisson fluctuations on scales larger than 
several arcmin (see Fig. 3). This means that sub-mm maps filtered with resolutions of several 
arcmin are a powerful tool to detect candidate proto-clusters of dusty galaxies. 

In other words, the relatively large beams used by CMB experiments collect photons 
from regions with Mpc physical sizes at z = 1-3 (at these redshifts the physical scale is 
~ 0.5 Mpc/arcmin) thus summing the contributions from tens of star-forming galaxies in over- 
dense regions. As first pointed out by [98] the galaxy clumps may then become detectable, 
even if individual galaxies are well below the detection limit. Unbiased searches of high -2 
protoclusters thus become possible, overcoming the need of targeting possible signposts of 
high density peaks (high -2 radiogalaxies or powerful QSOs), as done so far. 

It must be stressed, however, that over-densities are not necessarily proto-clusters. Some 
may be random alignments, projected on the plane of the sky, of large scale structures 
either unrelated or connected by a cosmic filament almost aligned with the line of sight. 
Quantitative predictions were worked out by [98]. Briefly, within the standard gravitational 
clustering scenario, we expect that the distribution of intensity peaks due to source over¬ 
densities has a very large variance, resulting from three contributions, discussed below [see 
102], The mean number of objects inside a volume V centered on a source is: 


(N)p = n f [l + £(r)] dV, 
Jv 


(4.1) 


where n is the mean source number density. The variance around (N) p is: 

((N-(N) p ) 2 ) p = (N) p + 

n 2 [ [ [((ri,r 2 ) + £(r-i 2 )-£(n)£(r 2 )] dV l dV 2 (4.2) 

JvJv 


with 

COr,r 2 ) = Q{z ) [£(ri)£(r 2 ) + £(ri)£(r 21 ) + £(ri 2 )£(r 2 )], (4.3) 

Q(z) being the amplitude of the three-point angular correlation function. 

For a survey with a Gaussian angular response function [eq. (2.8)] and £(r) = (ro/r) L8 
we have 

nJ 2 = n £(r) f{6) dV ~ 25.9 nr^ 8 [D A (z)9 h ]' /2 (4.4) 

Jv 
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where Da is the angular diameter distance and 


£(r,z) = 6 2 (-M ef f, z) (dm ( r j z) 


(4.5) 


wher b(M e g,z ) is the bias factor, M e s is the effective halo mass, £dm(u z) is the spatial 
correlation function of dark matter halos. The mean luminosity of the “clump” around the 
central source with luminosity L m is 


Ld(z) — L m (z)+ 

[dL'L f ^L',z) [ [1 + Z(r,z)\mdV 


'v 


(4.6) 


with variance: 


= <*i + / dL L &(L ,z) / [1 + £(r,z)} f(0)dV 
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where f(8) is the response function of the instrument 

m = e - ( e / 0 ) 2 / 2 

with 

FWHM 
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FWHM being the Full Width at Half Maximum of the instrument. 

The first contribution to the variance refers to the luminosity of the most luminous source, 
acting as a beacon signalling the presence of the proto-cluster; the second refers to the sum 
of luminosities of the surrounding sources, that sample in a different way the luminosity 
function; the third refers to the overdensity of neighbouring sources. 

All-sky surveys, such as the one planned by COrE+, are optimally suited to pick up 
the rare cases in which all contributions conspire to yield an exceptionally high luminosity 
of the clump. Particularly favourable situations may occur when the main source is strongly 
lensed or the clump is within a filament almost aligned with the line of sight, so that the 
contribution of neighbouring sources (seen in projection) to the observed signal can be very 
large. 

Under the assumption that the statistics of the matter density distribution can be 
described by a log-normal function [98], the probability distribution function of L c \ is: 


exp 
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The clump luminosity function T(L c i, 2 :) obeys the normalization condition: 


l (L d , z) L cl dL cl = f $(L, z) LcLL . (4.11) 

A search for such clumps has been carried out by [97] by investigating the nature of Planck 
ERCSC sources that lie within an area of ~ 90 deg 2 observed as part of the Herschel Multi¬ 
tiered Extragalactic Survey (HerMES). Four sources were found to be associated with over¬ 
densities of Herschel sources with sub-nun colours suggesting redshifts ~ 1-2. The estimated 
surface density of these candidate high -z clumps was found to be consistent with the predic¬ 
tions by [98], although with large uncertainties. 

However new calculations taking into account the more recent determinations of the 
redshift-dependent sub-millimeter luminosity functions (Negrello et al., in preparation) yield 
lower clump surface densities than estimated by [98], and fall short of accounting for the 
[97] result. A very similar conclusion was reached by [103] using a completely independent 
approach, coupling hydro-dynamical zoom-in simulations with the recently developed radia¬ 
tive transfer code GRASIL3D. The origin of the discrepancy is unclear. On one side it may 
suggest that models and simulations underestimate the IR luminosities of the clumps. Al¬ 
ternatively, the overdensities discovered by [97] might be not individual proto-clusters but 
positive fluctuations in the number of proto-clusters within the Planck beam. In both cases 
the study of these overdensities provides important clues on the early evolution of large scale 
structure. 


4.2 Power spectrum of the Cosmic Infrared Background (CIB) 

The better angular resolution, at high frequencies, of all COrE+ options, compared to Planck , 
will allow us to measure, in a uniform way, the CIB power spectrum over an unprecedented 
range of frequencies and of angular scales (from ~ 1 arcmin to tens of degrees), thus break¬ 
ing the degeneracy between the Poisson contribution and that of non-linear effects (one-halo 
term), that complicates the interpretation of Planck measurements, without resorting to ex¬ 
ternal ( Herschel ) data. Although there is good agreement between the latest determinations 
of the power spectra by Planck and Herschel [15, 16], at 545 GHz the Herschel data tend to 
be slightly but systematically higher, although compatible within the error bars. 

Accurate determinations of the CIB power spectrum at different frequencies provide 
on one side constraints on the evolution of the cosmic star formation density and, on the 
other side, on halo masses associated to sources of the CIB. Except on the smallest angular 
scales, the CIB power spectrum is determined by galaxy clustering. The currently standard 
approach to describe it starts from the consideration that galaxies are hosted by dark matter 
halos. The power spectrum of the galaxy distribution is parameterized as the sum of the 
one-halo term that dominates on small scales and depends on the distribution of galaxies 
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within the same halo, and the two-halo term that dominates on large scales and is related to 
correlations among different halos. 

The halo occupation distribution, which is a statistical description of how dark matter 
halos are populated with galaxies, is modelled using a central-satellite formalism [see, e.g., 
104]. The two-halo term is essentially proportional to the square of the effective bias factor 
times the linear theory power spectrum of dark matter halos. Since the bias factor is a steep 
function of the halo mass and redshift, the two-halo term is a sensitive probe of the effective 
halo mass of CIB sources (i.e. of galaxies that contribute most to the cosmic star formation) 
and of its evolution with redshift. 

The large number of COrE+ channels will make it possible to investigate in detail the 
decorrelation of power spectra measured in different frequency bands. A decorrelation is 
expected because the redshift distribution of sources of the CIB shifts to higher and higher 
redshifts with decreasing frequency [24, 69, 105, 106], and was observed [15, 16]. Accurate 
measurements of the CIB cross spectra for different frequency channels sets strong constraints 
on the frequency dependence of redshift distributions, hence on the evolution of the cosmic 
SFR. 

5 Polarization of extragalactic sources 
5.1 Radio sources 

Based on our calculations, we expect 90% completeness limits in polarized flux density that 
are substantially lower than the corresponding limits in total intensity, shown in Fig. 4, and 
not far from the 5cr noise levels. This is because, in the COrE+ frequency range, the sky 
is much less complex in polarization than in temperature: the free-free, the spinning dust, 
the SZ and the CIB emissions are either unpolarized or very weakly polarized; the CMB is 
also weakly polarized and, moreover, has a different polarization pattern (dominated by the 
E-mode) than foreground emissions, which have comparable E- and B-rnode polarization. 
The source confusion level is also low since, roughly, the rms polarization fluctuations, a p , 
are related to the fluctuations in total intensity, cr,, by a p ~ Bcr*, where II is the mean 
polarization degree [107]. [108] found, for flat-spectrum radio sources that dominate the 

counts at high radio frequencies, II ~ 0.026 at 18 GHz. 

Little is known on the polarization degree of dusty galaxies, but it is likely to be low 
because the complex structure of galactic magnetic fields with reversals along the line of sight 
and the disordered alignment of dust grains reduce the global polarized flux when integrated 
over the whole galaxy. The measurements at 850 /mr of M82 by [109] gave a global net 
polarization degree of only 0.4%. 

The estimated COrE+ 90% completeness limits in polarization for the 1 m option de¬ 
crease from ~ 25 mJy to ~ 6 nrJy as the frequency increases from 60 to 200 GHz and increase 
at higher frequencies up to ~ 23mJy at 350 GHz. For the 1.5 m option the 5 a noise levels in 
polarization vary from ~ llrnJy at 60 GHz to ~ 3mJy at 200 GHz, to ~ 8mJy at 350 GHz. 

Adopting the slope of 1.3 for the integral counts of radio sources in terms of polarized 
flux density (N(> S p ) oc S ~ L3 ) found by [108] at 18 GHz [see also 110], we find that COrE-f- 
will detect up to a factor of 40 (lm option) or of 160 (1.5 m option) more polarized radio 
sources than can be detected by Planck. In other words, COrE+ will provide blindly selected 
samples of thousands of polarized radio sources. We caution that these estimates are only 
tentative because of the complex spectral behaviour of the polarized flux density [108] that 
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make extrapolations in frequency quite uncertain. On the other hand, this uncertainty adds 
interest to COrE+ measurements. 

5.2 Dusty galaxies 

If the net polarization degree of M82 is typical of star forming galaxies, we expect the 
detection in polarization of ~ 7.5 sr -1 dusty galaxies for the lm option and of ~ 37sr^ 1 
for the 1.5 m option, corresponding to several tens or a few hundred galaxies, respectively, 
in the “extragalactic zone” (|6| > 30°). This will provide valuable constraints on the global 
properties of magnetic fields in star forming galaxies and on their relationships with SFRs. 


6 Discussion and conclusions 

COrE+ will perform substantially better than Planck also in the case of a telescope of a 
similar or somewhat smaller (1 m) size. There are two main reasons for this. COrE+ will 
be confusion limited up to the highest frequencies, in contrast to Planck. For example, at 
545 GHz (550/mi) the Planck beam has a FWHM = 3.8', while the diffraction limit is 1.5 / . 
The better resolution implies substantially deeper point source detection limits because the 
survey depth is mostly limited by fluctuations of sky signals. Realistic simulations give, at 
545 GHz, a 90% completeness level of 96mJy for the COrE+ 1.5 m option and of 141 mJy 
for the 1 m option. For comparison the corresponding flux density limit for the PCCS is of 
570 mJy. 

At low frequencies, where Planck resolution is also at the diffraction limit, the better 
performance by COrE+ is due to its lower instrumental noise. Source detection algorithms 
efficiently filter out diffuse emissions (CMB and Galaxy), whose power spectra sink down 
rapidly at the relevant resolutions, but suffer instrumental noise, which has a flat, roughly 
white noise, power spectrum. Our PSM-based simulations give, at 100 GHz, 90% complete¬ 
ness limits of 85 and 180mJy for the 1.5m and lm options, respectively, to be compared 
with the limit of 300 mJy for the PCCS. 

At higher frequencies, the Planck ERCSC has already provided large samples of dusty 
galaxies mostly at distances < 100 Mpc, thus offering the first opportunity to accurately 
determine the local luminosity function unaffected by cosmological evolution and to inves¬ 
tigate their properties. The PCCS reaches distances about a factor 1.55 larger, increasing 
the explored volume by about a factor of 3.7. COrE+ will reach still further, allowing us to 
explore a volume a factor of about 17 (1.5m option) or of about 8 (lm option) larger than 
the PCCS. Spectroscopic or photometric redshifts will be available for essentially all these 
galaxies thanks to the large area surveys from the ground and to the Euclid slitless spec¬ 
troscopy. This will make analyses of the properties of galaxies (luminosity functions, dust 
mass functions, star formation rate functions, dust temperature distributions, etc.) possible 
across the Hubble sequence. 

But the main improvement of COrE+ compared to Planck observations in total intensity 
will be in early galaxy evolution. The sub-mm surveys by COrE+ will cover the gap between 
the essentially Euclidean portion measured by Planck and the steep one measured by the 
deep SPT and Herschel surveys. The transition between the two portions provides a strong 
test for galaxy evolution models. Many of them predict it to be smoother than suggested by 
present data, but only the much better COrE+ statistics will make possible to reach firm 
conclusions. 
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Planck has already found high-redshift galaxies with extreme intrinsic IR luminosities, 
further boosted by strong gravitational lensing. Hyperluminous IR galaxies are very inter¬ 
esting in themselves since have always posed one of the strongest challenges to semi-analytic 
galaxy evolution models. Based on our simulations we expect that COrE+ will detect, at 
600 GHz, from about 321 sr -1 (lrn option) to about 992 sr" 1 (1.5 m option) strongly lensed 
galaxies, i.e. from a few to several thousands such objects in the high Galactic latitude 
sky. Such large samples are of extraordinary astrophysical value in many fields because the 
magnification makes galaxies that are intrinsically too faint become observable, while the 
corresponding increase of the apparent size makes it possible to measure the internal struc¬ 
ture to levels otherwise unattainable. Optical spectroscopy of galaxies acting as lenses can 
be exploited to measure the mass distribution of their dark matter halos as a function of 
redshift. This will allow a direct test on the evolution of large-scale structure. Samples of 
thousands of strongly lensed galaxies are essential for many other astrophysical and cosmo¬ 
logical applications, e.g. to measure cosmological parameters. Essentially all high -2 dusty 
galaxies will be strongly lensed. The other sources detected at high frequencies will be easily 
recognizable low -2 late-type galaxies, plus a small group of radio sources, also easily identi¬ 
fiable in low frequency radio catalogs. The selection of strongly lensed galaxies will then be 
extremely easy. 

At sub-mm wavelengths fluctuations on scales of several arcrnin are dominated by clus¬ 
tering of high-z massive galaxies, interpreted as progenitors of present day giant ellipticals. 
Sub-mm surveys are therefore optimally suited to look for proto-clusters at earlier redshifts 
than can be reached by optical/near-IR, X-ray, SZ surveys. Herschel observations of intensity 
peaks with red sub-mm colours, suggestive of high redshifts, in Planck high frequency maps 
have indeed shown that several of them represent the joint emission from multiple starbursts 
at high 2 , as expected for a young galaxy cluster. 

Again, COrE+ can do much better than Planck. Its high sensitivity sub-mm maps, 
filtered at several arcmin resolution corresponding to proto-clusters scales at 2 = 1-3, will 
be a powerful tool to carry out unbiased searches of the long-sought galaxy proto-clusters 
whose intergalactic gas has not yet reached the virial temperature (i.e. is not detectable in 
X-rays or via the SZ effect) and whose member galaxies are dust obscured, i.e. optically faint. 
Only very large area surveys, such as COrE+, can provide enough statistics on these rare 
objects. Higher sensitivity and higher angular resolution follow-up (sub-)mm measurements 
will assess the nature of candidate proto-clusters. This would help to solve a central problem 
in cosmology: how did the large scale structure of galaxies form? 

COrE+ will also measure, in a uniform way, the CIB power spectrum over an un¬ 
precedented range of frequencies and of angular scales (from ~ V to tens of degrees), thus 
breaking the degeneracy between the Poisson contribution and that of non-linear effects (one- 
halo term), that complicated the interpretation of Planck measurements, without resorting 
to external (e.g. Herschel or SPT) data. 

Moreover, COrE+ will extend the counts of radio sources compared to Planck , most 
notably at high frequencies, where very little is known. Above 217 GHz, the counts will be 
determined for the first time over a substantial flux density range with good statistics. This 
will enable the first investigation of the (sub)-mm spectral energy distribution (SED) and 
of the evolutionary properties of radio sources at (sub-)mm wavelengths. The vast majority 
of these sources are expected to be blazars, and the accurate determination of their spectra 
will allow us to understand how physical processes occurring along relativistic jets shape the 
SED. We will also get numerous samples of ‘extreme’ radio sources allowing us to investigate 
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the rich phenomenology of radio sources at (sub-)mm wavelengths. 

A spectacular advance will be made possible by the high COrE+ sensitivity in the blind 
detection of extragalactic sources in polarization. We estimate that COrE+ will detect up 
to a factor of 40 (lrn option) or of 160 (1.5 m option) more polarized radio sources than can 
be detected by Planck. We also expect, for the first time, blind samples of several tens (1 m 
option) to a few hundred (1.5 m option) of star forming galaxies detected in polarization in 
the “extragalactic zone”. These samples will open a new window for studies of the global 
properties of magnetic fields in star forming galaxies and of their relationships with SFRs. 
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